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ABSTRACT 


An  analytical  investigation  of  laminar  free  convection 
from  a  semi-infinite  vertical  plate  is  undertaken  with  the 
following  special  conditions  imposed.  The  plate  is  main¬ 
tained  at  a  constant  temperature  greater  than  that  of  the 
ambient  fluid.  At  an  arbitrary  distance  above  the  leading 
edge  a  constant  rate  of  blowing  or  suction  is  initiated 
extending  uniformly  along  the  remainder  of  the  plate.  The 
solution  was  obtained  using  the  method  of  matched  asymptotic 
expansions  whereby  an  inner  series  described  the  velocity 
and  temperature  profiles  near  the  discontinuity  in  the 
vicinity  of  surface  and  an  outer  series  approximated  the 
behavior  at  greater  distances  from  the  plate.  The  inner 
series  was  also  used  to  determine  the  heat  transfer  charac¬ 
teristics  at  the  surface  of  the  plate.  The  ambient  fluid 
was  assumed  to  be  air  with  a  Prandtl  number  of  0.72. 
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I.  INTRODUCTION 


The  study  of  free  convection  has  been  the  subject  of 
ongoing  research  over  the  past  century  both  from  an  analyti¬ 
cal  and  experimental  point  of  view.  An  excellent  overview 
of  the  progress  made  in  this  area  can  be  found  in  Ede  [1]. 

Two  outstanding  expositions  of  the  fundamentals  of  heat 
transfer  by  the  mechanism  of  free  convection  are  presented 
by  Gebhart  [2]  and  Moore  [3]. 

The  technique  of  solving  boundary  layer  problems  by 
employing  similarity  analysis  is  discussed  in  general  by 
Hansen  [4]  and  specific  applications  of  the  method  are  pre¬ 
sented  by  Yang  [5]. 

The  problem  discussed  in  this  paper  is  that  of  an  analyti¬ 
cal  investigation  of  the  effects  of  blowing  and  suction 
applied  at  a  specified  distance  above  the  leading  edge  of  a 
uniformly  heated  semi - inf inite  vertical  plate  on  the  heat 
transfer  characteristics  at  the  surface  of  the  plate. 

An  early  study  of  the  effects  of  mass  transfer  on  free 
convection  was  conducted  by  Eichhorn  [6]  who  showed  that 
similar  solutions  are  possible  for  blowing  or  suction  rate 
distributions  which  vary  as  the  distance  from  the  leading 
edge  of  the  plate  raised  to  a  power  related  to  the  variation 
in  surface  temperature. 


13 


1 


Since  the  prescribed  surface  temperature  itself  varied  as 
a  power  of  the  distance  from  the  leading  edge,  this  proved  to 
be  a  somewhat  restricted  case  of  limited  physical  application. 
The  case  of  a  uniformly  heated  plate  with  constant  blowing  or 
suction  along  the  entire  extent  of  the  surface  was  succinctly 
treated  by  Sparrow  and  Cess  [7]  who  extended  their  analysis 
to  include  independent  power  law  variations  of  both  surface 
temperature  and  mass  transfer.  Their  analysis  which  was  based 
on  an  asymptotic  series  expansion  of  the  stream  function  and 
temperature  distribution  evaluated  the  first  order  effects  of 
blowing  and  suction  on  free  convection.  An  extension  of  this 
work  to  include  possible  variations  in  the  fluid  properties 
was  conducted  by  Parikh,  Kays  and  Bershader  [8]  who  employed 
a  numerical  finite  difference  scheme  to  solve  the  resulting 
equations.  Their  numerical  predictions  were  in  close  agree¬ 
ment  with  interferometric  data  obtained  through  their  own 
experiments.  The  results  of  these  authors  and  especially 
those  reported  by  Sparrow  and  Cess  have  been  of  great  value 
in  affording  a  comparison  with  results  obtained  in  the 
present  study.  A  comprehensive  investigation  of  the  effects 
of  suction  and  blowing  on  heat  transfer  and  wall  skin  fric¬ 
tion  which  included  regions  on  the  plate  at  large  distances 
from  the  leading  edge  was  conducted  by  Merkin  [9].  In  that 
analysis  a  step-by-step  numerical  solution  of  the  boundary 
layer  equations  begins  at  the  leading  edge  and  proceeds 
systematically  up  the  plate  until  the  asymptotic  velocity 
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and  teraeperature  profiles  are  obtained  to  whatever  accuracy 
is  required. 

There  have  been  several  additional  studies  in  the  area 
of  free  convection  past  a  vertical  plate  with  various  condi¬ 
tions  of  wall  temperature  or  heat  flux  and  imposed  variations 
in  transpiration  at  the  surface  of  the  plate.  An  analysis 
which  transformed  the  laminar  boundary  layer  equations  re¬ 
sulting  from  arbitrary  distributions  in  both  wall  temperature 
and  blowing  into  ordinary  differential  equations  with  variable 
coefficients  and  constant  boundary  conditions  was  presented 
by  Vedhanayagan,  Altenkirch,  and  Eichhorn  [10]  .  In  the  area 
of  unsteady  free  convection  where  either  the  plate  surface 
temperature  or  transpiration  rate  varies  with  time  there  have 
been  at  least  four  notable  studies  during  the  past  twenty 
years.  Nanda  and  Sharma  [11]  undertook  analytical  inves¬ 
tigations  and  determined  that  similarity  solutions  were 
possible  for  the  following  two  cases:  surface  temperature 
varying  as  some  power  of  time  and  suction  velocity  varying 
as  time  to  the  minus  one-half  power;  and  the  second  case 
where  surface  temperature  is  an  exponential  or  periodic 
function  of  time  and  the  suction  velocity  is  constant.  The 
effects  of  mass  diffusion  on  the  unsteady  free  convective 
flow  past  a  simi-infinite  porous  plate  with  constant  suction 
were  studied  through  mathematical  analysis  by  Soundalgekar 
and  Wavre  [12]  .  With  the  assumption  that  the  plate  temperature 
oscillated  in  time  about  a  constant  mean,  approximate 
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solutions  were  derived  for  the  mean  flow,  transient  flow,  the 
amplitude  and  phase  of  the  skin  friction  and  the  rate  of  heat 
transfer.  The  effects  of  the  nondemensional  parameters  Grashof 
number,  Prandtl  number,  Eckert  number  and  Schmidt  number  were 
also  reported.  The  most  general  unsteady  situation  where  both 
wall  temperature  and  suction  velocity  are  periodic  functions 
of  time  was  considered  by  Lai  [13]  .  All  of  the  previous 
studies  cited  have  neglected  the  effects  of  viscous  dissipation; 
however,  this  effect  is  not  insignificant  when  the  natural 
convection  flow  field  is  very  large  or  the  flow  occurs  in  the 
presence  of  extremely  low  temperatures  or  very  high  gravity 
fields.  The  combined  effect  of  viscous  dissipation  with  the 
suction  velocity  and  plate  temperature  oscillating  at  the 
same  frequency  on  free  convection  past  a  vertical  plate  was 
analyzed  by  Soundalgekar  and  Pop  [14]  .  In  their  study  expres¬ 
sions  for  the  mean  velocity,  transient  velocity,  transient 
temperature  profiles,  amplitude  and  phase  of  skin  friction 
and  the  rate  of  heat  transfer  were  derived  and  the  results 
compared  to  those  reported  for  the  case  of  constant  suction 
velocity . 

A  review  of  the  literature  in  the  field  of  free  convection 
on  a  vertical  plate  clearly  indicates  a  desire  on  the  part 
of  researchers  to  investigate  perturbations  of  the  thermal 
and  momentum  boundary  layers  resulting  from  the  introduction 
of  temperature  or  velocity  discontinuities  along  the  surface 
of  the  plate  and  their  ultimate  effect  on  skin  friction  and 
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heat  transfer  at  the  wall.  The  specific  nature  of  the  dis¬ 
continuity  will  determine  the  manner  in  which  the  perturbation 
effects  propagate  through  the  boundary- layers  and;  therefore, 
will  often  indicate  a  preferred  method  of  analytical  solution. 
Although  there  are  several  possible  mathematical  techniques 
available  to  researchers  in  this  field,  two  fundamentally 
different  approaches  predeminate,  each  with  numerous  variations. 
An  example  of  the  finite  difference  method  of  solution  to  the 
boundary- layer  equations  resulting  from  a  discontinuity  in  the 
surface  wall  temperature  can  be  found  in  Na  [15] and  Meena  and 
Nath  [16]  .  This  technique  has  gained  in  acceptance  with  the 
advent  of  larger  and  faster  digital  computers.  The  other 
major  technique  involves  the  expansion  of  both  the  stream 
function  and  temperature  profile  in  an  asymptotic  series. 

The  original  proponent  of  this  method  was  Goldstein  who  pre¬ 
sented  the  details  of  the  asymptotic  series  technique, 
commonly  known  as  the  Goldstein  series,  in  [17],  After  re¬ 
ducing  the  set  of  governing  boundary- layer  equations  by  the 
introduction  of  a  stream  function  such  that: 


3i b 
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u  =  and  v 

two  new  variables  are  introduced  which  are  of  the  following 
general  form: 

1)  In  the  streamwise  direction:  £  =  x' 

2)  In  the  direction  normal  to  the  plate:  q  =  yx 
A  series  is  then  formed  in  powers  of  £  for  the  stream  function 
and  temperature  distribution  with  undetermined  coefficient 
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functions  f y ,  g^  respectively  which  depend  only  on  n.  The 
problem  is  then  finally  reduced  to  one  of  determining  the 
coefficient  functions;  the  solution  of  a  fifth  order  set  of 
linear  coupled  ordinary  differential  equations  for  each  pair 
of  coefficients  subject  to  the  appropriate  boundary  conditions. 
This  technique  is  especially  powerful  when  employed  to  inves¬ 
tigate  conditions  at  or  near  the  surface  of  the  plate.  In 
order  to  examine  the  laminar  free  convection  wake  above  a 
heated  vertical  plate  Yang  [18]  successfully  employed  this 
technique  to  obtain  detailed  velocity  and  temperature  profiles 
in  the  immediate  vicinity  of  the  trailing  edge  and  then 
expanded  these  results  to  the  rest  of  the  wake  by  an  integral 
solution.  Another  example  of  this  technique  is  found  in 
Kelleher  [19]  where  a  step  discontinuity  in  the  surface 
temperature  imposed  at  a  finite  distance  above  the  leading 
edge  of  the  plate  caused  a  perturbation  of  both  the  thermal 
and  momentum  boundary- layers .  An  inner  expansion  of  the 
Goldstein  type  was  employed  to  obtain  the  behavior  of  heat 
transfer  near  the  surface  of  the  plate  and  an  outer  asymptotic 
series  expansion  was  constructed  to  determine  the  velocity 
and  temperature  profiles  at  greater  distances  from  the  plate. 

A  similar  approach  is  employed  to  analytically  investigate 
the  present  problem  where  a  discontinuity  in  the  velocity 
boundary  condition  is  propagated  through  the  momentum  and 
thermal  boundary- layers .  Laminar  flow  is  assumed  as  are 
constant  transport  properties  of  the  ambient  fluid.  The 


primary  aim  of  the  study  is  the  determination  of  any  en¬ 
hancement  or  degradation  of  the  heat  transfer  rate  at  the 
surface  of  the  plate  compared  to  the  same  plate  without 
blowing  or  suction  over  the  entire  length. 
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II.  ANALYSIS 


A.  COORDINATE  SYSTEM 

The  problem  considered  in  this  paper  is  that  of  laminar 
free  convective  flow  from  a  uniformly  heated  semi- inf inite 
vertical  plate  with  leading  edge  oriented  downward  in  a 
standard  gravitational  field.  A  constant  rate  of  blowing  or 
suction  is  applied  at  the  surface  beginning  at  an  arbitrary 
distance,  L,  above  the  leading  edge.  A  rectangular  cartesian 
coordinate  system  with  origin  as  indicated  in  Figure  1  was 
employed  in  the  analysis.  As  is  indicated  throughout  the 
paper  various  transformations  of  the  basic  coordinates  are 
invoked  to  facilitate  the  solution  of  the  problem. 


B.  GOVERNING  EQUATIONS 

The  analysis  is  based  on  the  standard  Boussinesq  free 
convection  approximation  which  assumes  that  only  density  va¬ 
riations  giving  rise  to  buoyant  forces  are  considered  with 
the  effects  of  viscous  dissipation  neglected.  With  these 
standard  assumptions  the  two  dimensional  equations  of  con¬ 
tinuity,  momentum,  and  energy  boundary- layer  form  are: 


3u  f  3v 
3x  9y 


0 


(1) 


ulH+vl^=viiH+  gg(T  -  T J  (2) 

3x  3y  3y2 
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u  il  *  7  II  -  a 
3x  3y 


(3) 


3y 2 


where  T  is  the  constant  ambient  fluid  temperature.  These 
00 

equations  must  also  satisfy  the  following  physical  boundary 
conditions : 


At  y  =  0:  u  =  0 

7  3  7w 
v  =  0  , 

T  *  T 

w 


x>0 

- L<x<0 


(3a) 


and  as  y  -*■  00 : 


u  —  0 
T  =  T 


The  non- zero  normal  velocity  boundary  condition  at  x  =  0  and 
y  =  0  introduces  a  discontinuity  into  the  boundary- layer  equa¬ 
tions  which  has  a  primary  direct  effect  on  the  formation  of 
the  momentum  boundary- layer  and  a  secondary  influence  on  the 
development  of  the  thermal  boundary- layer .  A  positive  value 

for  v  will  represent  the  blowing  of  fluid  away  from  the  wall 
w 

and  a  negative  value  for  vw  will  represent  the  suction  of  fluid 
through  the  plate. 

In  order  to  facilitate  the  analysis,  the  variables  in 
equations  (l)-(3)  are  non-dimens ionali zed  according  to  the 
following  transformations: 
x 


x  = 


y-l 


(4a) 


(4b) 
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u 


(4c) 


V  = 


vL 

v 


(4d) 


9 


(4e) 


where  L  is  the  distance  above  the  leading  edge  where  suc¬ 
tion  or  blowing  commences,  and  is  the  constant  wall 
temperature.  The  equations  are  further  simplified  through 
the  introduction  of  the  conventional  stream  function  and 
by  stretching  y,  u,  and  v  as  follows: 


y*  - 


U 


*  s 


u(4Gr) 


-h  _  3</> 
3v* 


r*  * 


v(64Gr) 


-k  _  H 
3x 


(5a) 

(5b) 

(5c) 


where  the  Grashof  number,  Gr ,  is  evaluated  with  respect  to 
L,  the  length  scale  associated  with  the  discontinuity.  It 
should  be  noted  that  this  length  scale  is  not  the  conven¬ 
tional  one  used  in  free  convection  problems  where  distance 
from  the  leading  edge  is  normally  used  to  evaluate  Gr.  Now 
after  substituting  equations  (5)  into  equations  (l)-(3),  the 
continuity  equation  is  eliminated  and  the  momentum  and  energy 
equations  reduce  to  the  following: 


4 

3y* 


3x3y* 


„  _  3s*  .  n 

4  3*  3y*T  3y*T  +  9 


(6) 
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,  3»  39 

3y*  3x 


.  3i|>  30 
4  37  ly 


(7) 


]__  3  2  0 
Pr  3y*2 


The  non- zero  normal  velocity  component  imposed  at  the  surface 
of  the  plate  perturbs  the  governing  equations  (6)  and  (7) 
precluding  the  existence  of  a  similarity  solution.  There¬ 
fore,  an  approximate  solution  to  these  equations  will  be 
obtained  by  means  of  inner  and  outer  asymptotic  series  ex¬ 
pansions  for  both  i p  and  0.  The  inner  expansions  will  approxi¬ 
mate  local  behavior  near  the  surface  of  the  plate  in  the 
vicinity  of  the  discontinuity,  and  the  outer  expansions  will 
model  velocity  and  temperature  profiles  at  greater  distances 
from  the  plate. 

C.  BOUNDARY  AND  INITIAL  CONDITIONS 

The  nondimens ional  form  of  the  governing  equations  given 
in  (6)  and  (7)  must  satisfy  the  original  physical  boundary 
conditions  as  specified  in  (3a)  which  reduce  to  the  following 
conditions  in  terms  of  the  stream  function  ip  and  temperature 
distribution  0: 

At  y*  =  0:  u*  =  =  0 

V*  =  =  ±£*  ,  x>  0 

v*  =  'll  =  0,  - L<x<0  (8a) 

o  X 

0  =  1 

As  y*  ®:  u*  =  0 

0=0 
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where  e*  is  related  to  the  physical  blowing  or  suction  rate, 
v  ,  by  the  following: 

-v  L  , 

e*  =  (64Gr)  "**  (3b) 

Since  this  problem  involves  a  perturbation  of  both  the 
momentum  and  thermal  boundary- layers  which  began  developing 
at  the  leading  edge,  the  initial  conditions  at  the  origin  of 
the  discontinuity,  i.e.,  x  =  0,  must  match  these  free  con¬ 
vection  profiles.  For  a  semi- inf ini te  vertical  plate  main¬ 
tained  at  a  constant  temperature  the  following  ordinary 
differential  equations  describe  unperturbed  free  convection 


behavior : 


+  3FF"  -  2F'2  +  H  =  0 


H"  +  3FH '  =  0 


with  the  following  boundary  conditions: 

At  n=  0 :  F(0)  =  F' (0)  ;  H(0)  =  1 

As  p-*-00:  F'(°°)  =  0;  H(°°)  =  0 


where  F  is  the  similarity  stream  function  and  H  is  the 
dimensionless  temperature  distribution  and  primes  indicate 
derivatives  with  respect  to  the  free  convection  similarity 
variable  defined  as  follows: 


y  rGr* A 

Z  > 

(x+L)  4 


where  Gr*  is  evaluated  at  x+L  according  to  the  free  convection 
length  scale  which  has  its  origin  at  the  leading  edge.  It 
should  be  noted  that  all  references  or  comparisons  to  the 
unperturbed  free  convection  problem  employ  this  translated 
length  scale,  whereas  the  length  scale  for  the  perturbed 
problem  is  based  on  L  with  origin  at  x  =  0.  At  x  =0  equation 

(11)  reduces  to  p£  =  y* .  The  solution  to  the  equations  (9) 
are  well  known;  however,  an  excellent  summary  of  solutions 
for  this  and  other  more  general  cases  can  be  found  in  Yang 
[5]. 

Since  the  initial  length  of  the  plate  is  unperturbed,  the 
velocity  and  temperature  profiles  at  x  =  0  will  be  given  by 
the  following: 

u  =  (4Gr)^  F ’ (y*)  (12a) 

Q  =  H(y*)  (12b) 

Therefore,  the  initial  velocity  and  temperature  profiles 
originating  at  x  =  0  must  match  those  obtained  from  equations 

(12) .  The  details  of  expressing  F  and  H  in  a  form  which  will 
permit  this  matching  condition  to  be  mathematically  formulated 
are  presented  in  Appendix  A. 
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III.  METHOD  OF  SOLUTION 


A.  REDUCTION  OF  THE  GOVERNING  EQUATIONS 

The  nondimens ional  forms  of  the  governing  equations  for 
the  problem  under  consideration  are: 


Ty*  3x3yJ 


,  3^  324*  _  33^  . 

3x  3y*  3 


4  31»  li  .  .  3^  .  39  .  1  3 2  8  7 

Sy*  3x  H  3x  Sy*  “  Pr  3y^ 

In  order  to  apply  the  Goldstein  series  method  of  solution  to 
these  equations  suitable  expansions  of  ip  and  0  must  be  con¬ 
structed  of  the  following  form: 

1 

*  £*(^/<Pl)a  +  4>3Cc[Ck£kCn)]  (13) 

and  0  =  1  +  <J>4?d[£kgk(n)]  (14) 

where  $  =  <f>2xa;  n  =  <J>2y*x^;  the  s  are  possibly  functions  of 
e*;  and  the  a,  b,  c,  d  are  constants  to  be  determined.  The 
following  conditions  must  be  satisfied  which  impose  certain 
restrictions  on  the  <j>'s  and  exponent  constants: 

1)  The  u  velocities  from  the  free  convection  region 


(x<0)  and  from  the  perturbed  stream  function,  i.e., 

u*  =  Trr3F>  must  match  at  x  =  0. 
dy 

2)  The  temperature  profiles  must  match  at  x  =  0. 


3)  The  highest  order  derivative  of  the  coefficient  func¬ 
tions  must  appear  in  the  zeroth  order  expansions  which 


result  upon  substitution  of  equations  (13)  and  (14) 
into  equations  (6)  and  (7)  . 

4)  The  exponent  constants  in  the  expansions  referred  to 
in  3)  above  must  be  integers. 

The  details  of  applying  the  conditions  outlined  above  to 
determine  the  unknown  parameters  are  presented  in  Appendix  B 
along  with  the  boundary  conditions  which  must  be  satisfied 
by  the  coefficient  functions  f^  and  g^. 

Using  the  results  of  Appendix  B  equations  (13)  and  (14) 
are  transformed  into  the  following: 

*  =  +  C2Ukfk]  (15a) 

and  0  =  1  +  U£kgk]  (15b) 

where  £  *  (jx)1^  and  n  =  (|x)  " ly/ 3y*  =  y*/£  .  (16) 

The  momentum  and  energy  equations  (6)  and  (7)  are  transformed 
into : 

[(k  +  1)  Skfk]  Kkfkl  -  [(k  +  2)£kfk][£kf£] 

-  4e*£[5kf£]  -  [Ckfk' ]  +  5  +  C2[Ckgk]  (17) 

and  [(k  +  l)5kgk] [Ckf^]  -  [(k  +  2) £kfkJ  [£kg£] 

-  4e*C[ekg^l  =  ^pUkg£]  (18) 

where  the  repeated  index  k  indicates  the  standard  Einstein 
summation.  By  the  process  of  equating  coefficients  of  like 
powers  of  £  the  following  ordinary  differential  equations 
involving  the  coefficient  functions  are  obtained. 
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Momentum  Equation: 


Oth 

order : 

►+> 

+ 

IS) 

H} 

i 

f'2  =  0 

0  0  0 

• 

1st 

order : 

ft  t  .  +  2 f  f"  - 

1  8  1 

3f '  f '  +  3 f " f 

0  i  0  1 

=  -  4e  *  fj' 

-  1  (19) 

2nd 

order : 

f "  '  +  2f  f"  - 
2  0  2 

-  4e  *  f "  -  g 

1  0 

4f '  f  ’  +  4 f "  f 

0  2  0  2 

=  2  f  ' 2 
i 

-  3f  f" 

i  i 

Energy  Equation: 

Oth 

order : 

i-  ct"  +  2 f  g' 
Fr  °o  o° o 

-  f;s.  ■  ° 

1st 

order : 

i-  a"  +  2  f  g' 

Pr  ° i  si 

•  H'l,  *  3g;f, 

-  cr  f' 
»  1 

= -4e *g ' 

0 

(20) 

2nd 

order : 

1-  a"  +  2  f  g* 

Pr  s2  os2 

-  3ff2  +  4<j  *  f 
0S2  °0  2 

-  a  f  ' 
°0  2 

-  2g>f; 

-  3g'1fl  -  4e  *g 

t 

1  1 

From 

Appendix 

B  the  following  boundary  conditions 

must  be 

satisfied  by 

the  zeroth  order  momentum  and 

energy 

equations : 

At 

n  *  0  : 

f  CO)  -  f '  (0) 

0  0 

*  g0(0)  =  0 

As 

fl  ,  So 

rf0,  =  2a2  ’  “  “ 

■  bl 

where  a  =  *sF"(0)  and  b  *  H'(0). 

2  1 

The  solutions  of  the  zeroth  order  equations  are  now  straight¬ 
forward  and  yield  the  following: 

(21a) 
(21b) 

B.  DETERMINATION  OF  THE  HIGHER  ORDER  COEFFICIENT  FUNCTIONS 


fo  =  a2r'2 

and  go  =  b  ^  n 


The  first  and  higher  order  coefficient  functions  f^,  and 
for  the  stream  function  and  temperature  distribution  inner 
expansions  are  obtained  by  solving  the  sets  of  coupled  fifth 
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order  ordinary  differential  equations  (19)  and  (20)  .  These 
equations  were  obtained  by  equating  like  powers  of  £  in  the 
momentum  and  energy  expansion  defined  in  equations  (17)  and 
(18).  The  presence  of  e *  in  the  inhomogeneous  parts  of  the 
differential  equations  introduces  another  independent  para¬ 
meter  which  significantly  complicates  the  numerical  procedure 
employed  in  their  solution.  To  eliminate  this  problem  and 
isolate  the  effects  of  e*  the  coefficient  functions  can  be 


redefined  in  terms  of  universal  functions  which  depend  only 
on  the  Prandtl  number.  The  revised  forms  of  f^.  and  g ^  appear 
in  equations  (22),  presented  below. 

f  *  s*'Y 

*-i  =  z..  *  €*2.. 

f  «Y  +  e  *  Y  ♦  e*V  (22) 

1  1  20  *21  ‘22  y  J 

a  s  Z  +  £  *Z  +  £  *  2  Z 

°  2  2  0  21  22 


where  the  Y  ,  Z  ,  etc.  are  universal  functions  of  n  depend- 
10  10 

ing  only  on  the  Prandtl  number.  Since  the  first  and  higher 

order  equations  of  (19)  and  (20)  are  linear,  the  principle 

of  superposition  applies  and  the  equations  of  (22)  can  be 

substituted  into  (17)  and  (18)  and  separate  like  power  of 

e*  and  obtain  the  following  equations. 

The  first  order  momentum  equation  becomes: 

M  [Y  ]  =  -  1  (23a) 

1  i  o 

and  M  [Y  ]  *  -  4fM,  (23b) 

III  o 

where  M  is  the  linear  differential  operator  defined  by: 

i 

M  =  Y "  '  +  2f  Y"  -  3f  Y'  +  3f"Y  (23c) 

1  1  0  1  0  1  O'  1 
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The  first  order  energy  equation  becomes: 

E  [Z  ]  =  0 
1  i  o 

and  E  [Z  ]  =  -  4Z'  , 
xii  oo 

where  E  is  the  linear  differential  operator  defined  by: 
i 


E,  -  kz"  *  2f.z; 


2f '  Z  +  3Z '  Y  -  z  y; 
0  1  0  0  1  0  0  1 


The  second  order  momentum  equation  becomes: 

M  [Y  ]  =  2Y'2  -  3y  y"  -  Z  ; 

2  20  10  1010  00 

M  [Y  ]  *  4Y '  Y 1  -  3Y  Y" 

2  2  1  10  11  10  11 

and  M  [Y  ]  =  2Y ' 2  -  3Y  Y" 

222  11  1111  11 


3Y  Y"  -  4y" 
1110  10 


4Y' 


where  M  is  the  linear  differential  operator  defined  by: 
2 

M  =  Y"  *  +  2f  Y"  -  4f »Y'  +  4f’’Y  . 

2  2  02  02  02 

The  second  order  energy  equation  becomes : 

E  [Z  ]  =  2Y '  Z  -  3Y  Z'  ; 

2  2  0  1  0  1  0  1  0  1  0 


3Y  Z' 
10  11 


-  3Y  Z' 
11  10 


-  3Y  Z'  -  4Z' 


E  [Z  ]  -  2Y'  Z  +  2Y'  Z 
221  1011  1110 

-  4Z '  ; 

1  0 

and  E  [Z  ]  =  2Y '  Z 

222  1111  1111  11 

where  E  is  the  linear  differential  operator  defined  by: 
2 

E  =  4-  Z"  +  2f  Z*  -  3ffZ  +  4Z '  Y  -  Z  Y' 

2  2  02  02  22  002 

The  following  boundary  conditions  must  be  satisfied 
equations  (23) -(26): 

=  0  ; 


=  0 


(24a) 

(24b) 

(24c) 

(25a) 

(25b) 

(25c) 

( 2  5  d) 

(26a) 

(26b) 

(26c) 

(26d) 

by 


At 

ii 

o 

Y 

= 

Y  =  Y 

=  Y 

=  Y 

i 

0 

11  2  0 

2  1 

2  2 

Y  ' 

= 

Y'  =  y  ' 

=  Y' 

=  Y' 

l 

0 

11  2  0 

2  1 

2  2 

and 

Z 

= 

1-4 

II 

=  z 

=  Z 

i 

0 

11  2  0 

2  1 

2  2 

As 

n-**° : 

Y  i 

0 

-  i*  v  1 

0 

^ »  Y 

(27) 


i  l 
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and 


Y' 

2  0 

n3 
z 
n 

L 

2  0 

„3 


-H' (0) 


;  y' 

2  1 


0,  Y’ 

2  2 


0; 


1  0 
2 


0,  z 


2  1 


0; 


(28) 


0,  Z  •>  0,  Z 
2  1  2  2 


which  follow  directly  from  the  coefficient  function  boundary 
conditions  given  in  equations  (B-33),  (B-35)  and  B -  36 )  . 

Although  they  were  not  used  in  the  present  analysis,  the 
third  order  universal  function  equations  are  presented  in 
Appendix  C. 


C.  CONSTRUCTION  OF  THE  OUTER  SERIES  EXPANSIONS 

The  inner  series  developed  in  the  preceding  section  is 
useful  only  in  approximating  local  behavior  at  the  surface  in 
the  vicinity  of  the  discontinuity.  In  order  to  describe  the 
temperature  and  velocity  profiles  at  greater  distances  from 
the  plate  another  series  is  required  which  will  satisfy  the 
asymptotic  boundary  conditions  and  also  match  the  inner  series 
at  some  intermediate  distance  from  the  surface.  According 
to  the  method  of  Goldstein  [17]  the  outer  series  expansions 
for  the  stream  function  and  temperature  distribution  should 
be  of  the  following  form: 

-i|>  (Y*)  +  (Y*)  +  ~  Zzl>  (Y*)  +  lT  (Y*)  (29) 

and 

0  -  0  (Y*)  +  ce  (Y*)  +  jr  z2q  (Y*)  ♦  We  (Y*)  (30) 

0  1  2  •  2  J :  3 
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Since  as  C^O  the  original  initial  conditions  must  be  satisfied, 
the  following  equalities  must  hold;  thus  defining  the  zeroth 
order  coefficients: 


P  (.v*)  =  F(y*) 
o 


9  Cy*)  =  H  (y*) 
o 

where  F  and  H  are  the  free  convection  similarity  functions  for 
the  stream  function  and  temperature  distribution,  respectively. 
This  condition  is  merely  a  mathematical  restatement  of  the  phy¬ 
sical  condition  which  requires  that  the  temperature  and  velocity 
profiles  at  the  discontinuity,  C=  0,  match  those  generated  by 
the  free  convection  process  in  the  region  between  the  leading 
edge  and  the  discontinuity. 

Now  the  asymptotic  outer  series  froms  (29)  and  (30)  are 
substituted  into  the  governing  equations  (6)  and  (7) .  Once  again 
through  the  process  of  collecting  terms  which  are  coefficients 
of  like  powers  of  £,  the  following  set  of  linear  ordinary  dif¬ 
ferential  equations  is  obtained  from  the  momentum  equation: 


p  •  P  •  -  p  "  p  -  o 
0  1  0  1 


p  »  p  '  +  p  '  2  —  p  •  p”  -  p  p  "  =  0 
0  2  1  11  2  2 


(32a) 


(32b) 


1  4>  'ip  *  +  p  'p  '  *  i-  p  'p  ’  -  4  p  p  ”  -  p  p  ,?  -i  p  p  ”  =  9  +p  "  ’  (32  c) 

2  0  3  12  221  212  2  1  2  3  0  00 

and  from  the  energy  equation  the  following  set  of  algebraic 


equations 


p  '  9  -  p  9  '  =0 

0  1  10 

p  <  e  +  p  *  0  -  p  0  •  -  p  0  ■  =  0 

0  2  11  11  2  0 


(33a) 

(33b) 
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i'Ke+ifj'e  +  yire  -  i  *  9'  -  p  e 1  -  i  *  e '  =  i-  e  "  (33c) 


0  3  12 


The  details  of  the  solution  of  equations  (32)  and  (33)  can  be 
found  in  Appendix  D  and  the  results  are  summarized  below: 


*2  - 

2 

p  =  6\p 


2y*'P  '  +  -n-  P ' 

0  Za  o 


9=0 

l 


C 

1  Q  I 

Ta~  eo 

2 


9  ^  (jj-  -  2y*)  9  ' 

3  ^a.  o 

2 

With  these  coefficient  functions  determined,  velocity  and 
temperature  profiles  for  large  y*  can  be  generated  according 
to  the  following  expansions: 

C 

u*  =  fy*  3  F’(y*)  +  [4i-F"Cy*)]C2  +  [|f'(v*) 

,  D  2  (36) 

+  s-(2i~  *  2y*)F"(y*)K3 


C  ,  D 

9  =  H  (y*)  +  4~«'  (y*)  C2  +  '  2y*)H’(y*)53  (37) 

2  2 

D.  NUMERICAL  SOLUTION  OF  UNIVERSAL  EQUATIONS 

The  universal  equations  (23) -(26)  with  boundary  conditions 
(27)  and  (28)  were  solved  numerically  with  an  IBM  3033  series 
digital  computer.  Double  precision  machine  arithmetic  was 
used  for  all  computations.  A  least  squares  correction  method 


developed  by  Nachsteim  and  Swigert  {20]  was  used  as  shooting 
point  routine  to  determine  the  unknown  starting  values  for 
each  set  of  fifth  order  two  point  boundary  value  equations. 

A  standard  fourth  order  Runge-Kutta  initial  value  integration 
routine  was  employed  to  tabulate  values  for  the  free  convec¬ 
tion  functions  F,  F',  F",  H  and  H ' ;  the  zeroth  order  inner 
series  coefficient  functions  f  ,  f',  f",  g  and  g';  and  the 

O  0  O  O  0 

higher  order  universal  functions  Y  ,  Y  ,  Y  ,  etc.  These 

10  11  2  0 

tabulated  values  were  then  used  to  generate  the  inner  and 
outer  velocity  and  temperature  profiles  as  well  as  the  local 
Nusselt  number  at  the  surface  of  the  plate. 
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IV.  RESULTS 


A.  PRESENTATION  OF  RESULTS 

1 .  Velocity  and  Temperature  Profiles 
a.  Inner  Series  Profiles 

The  inner  series  expansions  for  velocity  and 
temperature  were  computed  form  the  tabulated  values  of  the 
universal  functions  which  were  obtained  form  the  numerical 
solution  of  the  fifth  order  ordinary  differential  equations 
defined  by  equations  (23)  through  (25).  The  velocity  profiles 
were  obtained  for  specified  values  of  e*  the  blowing  or  suction 
parameter  according  to  the  following  equation: 

u*  =  £[f  +  (Y '  +  e  *Y '  K  +  (Y*  +  e*Y  ’  +  e*2Y'  K2]  (38) 

0  10  11  20  21  22 

where  the  f ' ,  Y'  ,  Y'  ,  etc.,  are  functions  of  n.  In  order 
0  10  11 

to  facilitate  comparison  and  matching  of  the  inner  and  outer 
profiles  the  nondimens ional  velocity  u*  was  plotted  as  a  func¬ 
tion  of  y*  where  the  relationship  between  y*  and  n  is:  y*=5*n. 
Numerous  sets  of  profiles  were  generated  for  analysis  and 
comparison;  however,  only  the  profiles  corresponding  to  the 
following  parameters  are  presented  in  this  paper: 

1)  Blowing  and  suction  parameter,  e*,  of  magni¬ 
tudes  ±0.05  and  ±0.25  where  positive  values  cor¬ 
respond  to  suction  and  negative  values  to  blowing. 

2)  Streamwise  nondimensional  parameter,  £,  corres¬ 
ponding  to  x/L  ratios  of  .25,  .50,  and  1.0. 


The  plots  of  u*  versus  y*  for  these  parameters  are  presented 
in  Figures  2-13  and  for  ease  of  comparison  both  inner  and  outer 
profiles  are  included  on  the  same  plot. 


The  temperature  profiles  generated  by  the  inner 

series  expansions  were  obtained  from  the  following  equation: 

9  -  1  +  5[g  +  (3  +  +  (Z  +  e*Z  +  e*2Z  K2] 

0  10  11  20  21  22 

(39) 


where  the  g  ,  Z  ,  Z  ,  etc  are  functions  of  n;  however,  as 
0  10  11 

was  the  case  with  the  velocity  profiles,  9  was  plotted  as  a 
function  of  y*  with  y*  =  £*n*  The  same  parameters  were  used 
in  plotting  6  versus  y*  and  the  outer  profiles  were  included 
for  comparison.  These  plots  appear  in  Figures  14-25. 
b.  Outer  Series  Profiles 

The  outer  series  expansions  which  were  constructed 
in  the  preceding  chapter  approximated  the  behavior  of  the  ve¬ 
locity  and  temperature  distributions  at  greater  distances  from 
the  plate.  The  outer  velocity  profiles  are  given  by  equation 
(36)  which  is  written  below  for  convenience: 


C  , 

u*  =  F '  (y  *)  ♦  [jL-  F"  (y  *)  ]  £2  + 

[|F '  <y  *)  ♦  g-  (f|  -  Zy  *)  F”  (y  *)  ]  C  3  (36) 

where  F'  and  F"  are  the  first  and  second  derivatives  of  the 
free  convection  stream  function  evaluated  at  the  location  of 
the  discontinuity,  x  =  0.  They  were  numerically  calculated 
for  values  of  the  free  convection  similarity  variable  which 


is  identically  equal  to  y*  at  x  =  0.  The  constant  a  =*F"(0), 

2  - 

and  the  constants  C  and  D  are  computed  numerically  from 

12 

equations  (D-16)  and  (D-17)  respectively .  The  outer  velocity 
profiles  corresponding  to  the  same  e*  and  £  parameters  as  were 
specified  for  the  inner  profiles  are  shown  in  Figures  2  through 
13. 

The  outer  expansion  which  defines  the  temperature 

profile  valid  for  large  values  of  y*  was  given  in  equation  (37) 

and  is  repeated  below: 

C  ,  D 

e  -  H(y*)  +  ji—  H'*y*U2  +  i(-2-|-  -  2y*)H’(y*K3  (37) 

2  2 

where  the  constants  C  ,  D  ,  and  a  ,  are  defined  as  before,  and, 

1  2  2 

the  H  and  H'  are  the  free  convection  temperature  distribution 
and  its  derivative  evaluated  at  x  =  0  where  y*  is  the  free 
convection  similarity  variable.  Once  again  the  plots  of  9 
versus  y*  for  the  same  parameters  as  for  the  inner  profiles, 
are  presented  in  Figures  14  through  25. 

2 .  Heat  Transfer  Characteristics  at  the  Surface 

The  standard  method  of  characterising  the  heat  transfer 
at  the  surface  of  a  heated  plate  is  to  calculate  the  local 
Nusselt  number  at  the  point  of  interest.  In  this  particular 
problem  the  Nusselt  number  at  any  location  above  the  disconti¬ 
nuity,  x>0,  is  given  by  the  following: 

Nu  =  ± -  ^I|_  (40) 

(tw-t»)  y=  0 
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the  value  of  n  chosen  as  the  end-point.  A  summary  of  the 
effects  of  varying  the  end-point  on  both  the  starting  values 
and  outer  series  constants  is  presented  in  Table  I.  For  each 
of  the  end-point  values,  nencj»  listed  in  the  table,  a  complete 
set  of  tabular  profiles  were  generated  and  qualitatively  com¬ 
pared  to  determine  the  degree  of  matching  evidenced.  The 
results  of  this  process  revealed  that  the  smaller  the  value  of 
nend’  t*ie  better  was  the  overall  matching  of  the  profiles. 
However,  for  Fenj  <2.0  no  further  improvement  in  the  matching 
was  observed.  Also,  the  boundary  conditions  at  n  =  00  were 
satisfied  to  equal  degrees  of  accuracy  for  each  nen(j  selected. 
Therefore,  all  of  the  numerical  results  presented  were  based 
on  an  end-point  value  of  nenj  =  2.0.  It  was  also  observed 
that  the  outer  profiles  were  quite  insensitive  to  variations 

in  the  constants  C  ,  C  ,  etc,;  and  that  for  a  fixed  e*  and 

10  11 

£  the  parameter  Nu(Gr/4) "  1  was  insensitive  to  change:  in 
the  initial  values  which  were  caused  by  changing  nend. 
b.  The  Inner-Outer  Profile  Match  Point 

In  addition  to  the  qualitative  numerical  procedure 
outlined  above  to  effect  a  matching  of  the  inner  and  outer 
profiles,  another  approach  was  undertaken  based  on  the  appli¬ 
cation  of  an  overall  energy  balance  to  the  problem  under  con¬ 
sideration.  The  principle  of  conservation  of  energy  requires 
that  all  energy  entering  a  specified  control  volume  must  leave 
that  volume  providing  there  is  no  storage  of  energy  within  the 
volume.  In  the  case  of  the  present  problem  all  of  the  energy 


TABLE  I 


SUMMARY  OF  STARTING  VALUES  AND  MATCHING  CONSTANTS 


nend 

2.0 

3.0 

5.0 

7.0 

10.0 

Y"  (0) 

1  0 

IO*9 

10'9 

10"9 

0 

0 

ZM  (0) 

1  o  v  ' 

10'9 

10'9 

IO*9 

0 

0 

Y"  (0) 
i  l  v  ' 

2.5787 

2.8611 

3.0893 

3.1745 

3.2263 

Z*  (0) 

lr  1 

-1.5531 

-1.7623 

-1.9329 

-1.9994 

-2.0403 

Y"  (0) 
2  0^ 

io"4 

10"  4 

io"4 

io"5 

io" 1  ° 

Z'  (0) 

io"4 

io" 4 

io"5 

io"5 

io"8 

Y”  (0) 

2P  1 

-2.5437 

-3.3282 

-4.0277 

-4.3080 

-4.4860 

Z'  (0) 

.12905 

. 21959 

.345070 

.  4034 

.44235 

Y"  (0) 

2.15407 

2.3655 

2.7645 

2.9743 

3.1229 

Z'  (01 
2  2  V 

-.96015 

-.9905 

-1.18464 

-1.3222 

-1.4087 

C 

1  0 

0 

0 

0 

0 

0 

C 

i  i 

1.605 

1.636 

1.622 

1.615 

1.61 

D 

2  0 

0 

0 

0 

0 

0 

D 

2  1 

3.2534 

4.403 

7.022 

9.875 

14.289 

D 

2  2 

-2.8477 

-2.9755 

-2.9811 

-2.8945 

-2.788 
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entering  the  system  at  the  surface  of  the  plate  must  exit  in 
an  upward  direction  parallel  to  the  plate.  The  mathematical 
formulation  of  this  requirement  is  presented  below  in  integral 
form. 

f  cp (T  -  T  )udy  =  r  q"(x)dx  (42) 

-'o  P  vi  q  0 

The  left  hand  integral  must  be  divided  into  two  separate  in¬ 
tegrals  since  there  exists  an  inner  and  outer  region  where 
the  velocity  and  temperature  behavior  are  described  by  differ¬ 
ent  expans  ions -- the  inner  and  outer  profiles.  The  dividing 
line  between  these  two  regions  must,  therefore,  define  the 
match  point  of  the  inner  and  outer  profiles.  The  details  of 
evaluating  the  integral  equation  (42)  are  presented  in 
Appendix  E.  The  only  quantative  result  of  the  energy  balance 
analysis  was  that  an  exact  balance  could  be  achieved  only  if 
the  matching  point  between  the  inner  and  outer  expansions  is 
taken  to  be  y*  =  0.  That  is  the  integral  in  (42)  must  be 
evaluated  using  the  outer  series  to  account  for  energy  leaving 
the  volume — the  left  side- -and  the  inner  series  to  account 
for  the  thermal  energy  entering  the  volume  at  the  surface  of 
the  plate.  As  the  profiles  themselves  suggest,  the  outer 
profiles  when  viewed  in  an  integrated  or  global  sense  more 
nearly  satisfy  the  energy  balance,  although  in  a  pointwise  sense 
they  clearly  do  not,  especially  at  the  surface  of  the  plate. 
Since  behavior  at  the  surface  of  the  plate  whether  with  respect 
to  heat  transfer  or  shear  stress  is  of  primary  interest,  from 
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an  engineering  standpoint,  knowledge  of  the  exact  location  of 
the  match  point  is  not  essential,  since  the  inner  series  ex¬ 
pansions  are  always  used  to  furnish  results  at  the  surface. 

2 .  Velocity  Profiles 

Before  analyzing  the  velocity  profiles  some  method  of 

characterizing  the  relative  magnitude  of  the  blowing  or  suction 

velocity  should  be  established.  To  accomplish  this  v^  will  be 

compared  to  the  maximum  value  of  the  u-component  of  velocity 

in  the  case  of  unperturbed  free  convection  evaluated  at  the 

location  of  the  discontinuity,  x  =  0,  which  will  be  designated 

by  umax-  The  ratio  vw/umax  which  characterises  this  comparison 

is  reduced  to  the  following  nondimensional  form  through 

appropriate  substitutions: 
v  , 

^ -  =  5.124  e  *Gr  (43) 

umax 

For  a  Gr  s  10®  for  which  the  laminar  flow  regime  is  still 
applicable  the  relationship  expressed  in  (43)  becomes: 

v 

—  =  .0512  e*  (44) 

u 

max 


From  this  correlation  the  blowing  and  suction  parameters  used 
in  the  present  problem  correspond  to  the  following  specific 
relationships : 

e*  =0.05  implies  v  =  .00256  u  , 
r  w  max’ 

and  _  _ 

e*  =0.25  implies  v  =  .0128  u 

^  w  max 
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These  values  for  v^  are  consistent  with  the  original  assumption 
underlying  the  governing  boundary  layer  equations. 

A  careful  review  of  the  inner  profiles  plotted  in 
Figures  2-13  reveals  the  following  features  about  velocity  be¬ 
havior  with  suction  or  blowing: 

1)  As  suction  is  increased  at  a  specified  value 
of  the  location  of  the  maximum  u*  is  drawn 
in  closer  to  the  plate  and  its  magnitude  is 
decreased . 

2)  At  a  given  value  of  suction  the  location  of 
the  maximum  u*  is  moved  closer  to  the  plate 
and  the  magnitude  of  u*  is  decreased  as  the 
value  of  ^  is  increased. 

3)  The  results  for  blowing  exhibit  the  reverse 
effects  on  the  location  and  magnitude  of  u* 
as  those  described  in  1)  and  2)  for  suction. 

Numerical  data  on  the  effects  of  suction  and  blowing  on 
the  location  and  magnitude  of  the  maximum  u*  for  the  inner 
velocity  profiles  as  well  the  corresponding  data  for  the 
unperturbed  free  convection  case  is  presented  in  Table  II. 

Comparison  of  the  suction  data  with  the  free  convection 
data  reveals  that  the  location  of  the  maximum  u*  is  shifted 


closer  to  the  plate  and  also  that  its  magnitude  is  reduced  by 
15%  when  £  =0.25  and  by  11%  when  £  =  1.0  with  a  suction  mag¬ 
nitude  of  e *  =0.05.  When  the  suction  is  increased  to  e*=0.25 


the  effects  are  correspondingly  amplified.  Although  this  is 
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TABLE  II 


LOCATION  OF  INNER  PROFILE  MAXIMUM  VELOCITIES 


'x 


L  .25 

.50 

1.0 

Free  Convection: 

y* 

1.020 

1.070 

u* 

max 

.  3088 

.  3383 

.  3907 

Suction  e  *  =  0 .05 

y* 

.  8356 

.  8221 

.7314 

u* 

max 

.  2623 

.2563 

.2444 

Suction  e*=0.25 

y* 

.  6868 

.6634 

.4906 

u* 

max 

.2426 

.2194 

.1862 

Blowing  £*=-0.05 

y* 

.  8929 

.9086 

.9267 

u* 

max 

.2665 

.  2717 

.2851 

.9780  .9951  1.036 

,2589  .2745  .3333 


Blowing  e*=-0.25 

y* 

Uraax 


the  behavior  which  would  be  expected  for  suction,  the  percent 
reduction  in  maximum  velocity  predicted  by  the  inner  expansion 
is  probably  too  great  considering  the  relative  magnitude  of 
suction  applied.  In  the  case  of  blowing  the  inner  series 
fails  to  produce  results  consistent  with  expected  physical 
behavior  since  neither  the  location  of  the  maximum  u*  is  moved 
farther  from  the  plate  nor  is  the  magnitude  of  u*  increased 
above  that  observed  for  free  convection.  Also  the  degree  of 
matching  between  the  inner  and  outer  profiles  lessens  as  the 
suction  magnitude  or  distance  from  the  discontinuity  is 
increased . 

As  indicated  above  the  results  obtained  in  the  case 
of  suction  are  consistent  with  the  physical  situation  since 
in  this  case  heated  fluid  is  drawn  closer  to  the  plate  where 
viscous  forces  tend  to  reduce  the  maximum  velocity  in  the 
streamwise  direction.  In  the  case  of  blowing  heated  fluid 
is  pushed  farther  from  the  surface  away  from  the  influence 
of  viscous  forces  and  into  the  region  where  buoyancy  forces 
act  to  accelerate  the  flow  in  the  streamwise  direction  there¬ 
by  increasing  the  maximum  u*.  In  addition  for  a  specific 
value  of  either  suction  or  blowing  the  effects  described 
above  will  be  accentuated  as  the  distance  from  the  disconti- 
nuity  is  increased;  i.e.;  for  increasing  values  of  j-  or  £, 
since  the  blowing  or  suction  effects  will  have  had  a  greater 
opportunity  to  predominate  over  unperturbed  free  convection. 
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The  lack  of  positive  correlation  in  the  case  of  blowing  versus 
unperturbed  free  convection  might  be  attributed  to  any  of  the 
following : 

1)  Insufficient  terms  in  the  inner  series  expan¬ 
sion;  since  the  signs  of  the  coefficient 
functions  alternate  in  a  pairwise  manner, 
addition  of  another  term  would  possible  serve 
to  effect  a  better  correlation. 

2)  The  effects  of  blowing  become  significant  in 
the  region  where  the  inner  series  no  longer 
provides  a  valid  representation  of  the  velocity 
behavior . 

The  effects  of  suction  and  blowing  on  the  outer  velocity 
profiles  are  graphically  presented  in  Figures  2-13;  however, 
the  specific  effects  on  the  location  and  magnitude  of  the 
maximum  value  of  u*  are  given  in  Table  III  using  the  same 
format  as  Table  II  for  the  inner  profiles.  The  results  in 
the  case  of  suction  are  somewhat  inconsistent  with  the  behavior 
expected  relative  to  the  free  convection  case.  The  location 
of  the  maximum  velocity  does  not  move  closer  to  the  plate, 
and  its  magnitude  does  not  decrease  when  compared  to  the  free 
convection  profiles.  This  anomaly  suggests  that  the  effects 
of  suction  on  velocity  are  more  significant  in  the  region 
close  to  the  surface  of  the  plate  where  the  behavior  is  better 
approximated  by  the  inner  series  expansion.  Conversely,  in 
the  blowing  case  the  results  correlate  well  with  expected 
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TABLE  III 


LOCATION  OF  OUTER  PROFILE  MAXIMUM  VELOCITIES 


x. 

L'  .25 

.50 

1.0 

Free  Convection: 

y* 

1.020 

1.070 

1.150 

umax 

.  3088 

.  3383 

.  3907 

Suction  e*=0.05 

y* 

1.053 

1.1588 

1.3628 

u* 

max 

.  3110 

.3466 

.4212 

Suction  e*=0.25 

y* 

.9386 

.9519 

1.1072 

Umax 

.  3117 

.  3468 

.4160 

Blowing  e*=-0.05 

y* 

1.110 

1.240 

1.472 

u* 

max 

.  312 

.  350 

.433 

Blowing  £*=-0.25 

y* 

1.206 

1.370 

1.617 

u* 

max 

.  317 

.364 

.466 
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behavior.  Relative  to  the  free  convection  case  the  maximum 
value  of  u*  moves  away  from  the  plate,  and  its  magnitude  in¬ 
creases  from  1.0%  to  10.8%,  with  a  blowing  parameter  e*=-.05, 

as  —  increases  from  .25  to  1.0.  When  the  blowing  magnitude 

L. 

is  e*  =  -0.25,  the  value  of  u*  increases  from  2.7%  to  19.21 

max 

as  increases  from  .25  to  1.0.  Whether  or  not  this  increase 
is  excessive  must  await  comparison  with  experimentally  obtained 
data.  In  any  case  it  appears  that  a  very  small  perturbation 
of  the  momentum  boundary  layer  at  the  wall,  produces  relatively 
significant  changes  at  increased  distances  from  the  plate. 

3.  Temperature  Profiles 

The  effects  of  blowing  and  suction  on  the  temperature 
profiles  are  not  as  visible  in  the  plots  Figures  14-25,  as 
was  the  case  for  velocity  since  there  are  no  maximum  velocity 
peaks  to  observe.  However,  the  same  general  pattern  emerges 
as  did  with  the  velocity  profiles,  and  key  points  are  summarized 
below  for  the  inner  temperature  profiles: 

1)  For  a  given  location  above  the  discontinuity  the 
inner  profile  pivots  clockwise  about  0  =  1.0  as 
the  magnitude  of  suction  is  increased,  this  re¬ 
presents  an  increase  in  the  gradient  which  is 
characteristic  of  an  increased  heat  transfer  rate. 

2)  For  giv-'.n  magnitude  of  suction  there  is  a  marked 
increase  in  gradient  as  the  distance  from  the 
discontinuity  in  the  streamwise  direction  is  in¬ 
creased,  i.e.,  the  value  of  £  increases. 
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3)  Relative  to  the  unperturbed  free  convection  profiles 

Y 

the  effect  of  suction  for  a  given  value  of  j-  is  to 
increase  the  gradient. 

4)  For  blowing  the  reverse  effects  of  those  outlined 
in  1)  through  3)  above  are  observed  with  the  addi¬ 
tional  feature  that  the  inner  profiles  are  bowed 
out  as  well  decreasing  their  gradient.  This  be¬ 
havior  characterizes  a  decreased  heat  transfer 
rate . 

These  observed  phenomenon  correspond  with  what  physical 
intuition  suggests  should  occur.  In  the  case  of  suction  heated 
fluid  is  moved  closer  to  the  wall  thereby  reducing  the  thick¬ 
ness  of  the  thermal  boundary - layer  and  enhancing  the  rate  of 
heat  transfer  as  is  evidenced  by  the  increase  in  gradient.  The 
opposite  situation  in  which  blowing  forces  heated  fluid 
farther  from  the  plate  causes  the  thermal  boundary- layer  to 
thicken.  This  creates  an  insulating  effect  which  serves  to 
reduce  the  rate  of  heat  transfer  evidenced  by  a  decrease  in 
gradient . 

The  outer  temperature  profiles  again  reflect  the 
distinctive  shape  of  unperturbed  free  convection  profiles. 

In  the  case  of  suction  their  behavior  is  inconsistent  and 
they  do  not  behave  relative  to  the  free  convection  profiles 
as  do  the  inner  profiles.  This  again  suggests  that  the  effects 
of  suction  are  most  evident  near  the  surface  of  the  plate 
where  the  inner  series  is  a  closer  approximation  of  the  actual 
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behavior.  For  the  case  of  blowing  the  outer  profiles  consis¬ 
tently  move  away  from  the  plate  with  increases  in  either 
magnitude  or  distance  above  the  discontinuity.  Also,  rela¬ 
tive  to  the  free  convection  profiles,  they  behave  as  expected 
moving  farther  out  from  the  surface  as  the  magnitude  of 
blowing  is  increased  for  a  specified  value  of  This  suggests 
that  the  outer  profiles  are  indeed  a  good  approximation  of 
temperature  behavior  at  large  distances  from  the  plate. 

4 .  Heat  Transfer  at  the  Surface  of  the  Plate 

In  order  to  analyze  the  effects  of  discontinuous  suc¬ 
tion  and  blowing  on  the  rate  of  heat  transfer  at  the  plate  a 
graphical  comparison  was  made  of  the  following  three  physical 
situations : 

1)  Unperturbed  free  convection  from  an  isothermal 
plate . 

2)  Free  convection  from  an  isothermal  plate  with 
continuous  suction  or  blowing  applied  over  the 
entire  length,  (from  the  study  by  Sparrow  and 
Cess  [7]  . 

3)  Free  convection  from  an  isothermal  plate  with 
discontinuous  suction  or  blowing  applied  at  a 
distance  L  above  the  leading  edge. 

For  each  of  these  cases  the  heat  transfer  data  was  correlated 
according  to  following  nondimensional  parameter: 

Nu  (Gr/4) "  h  (46) 
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I 


where  the  Nusselt  number  in  each  case  Nu  ,  Nu,  and  Nu  was 

o  1 

evaluated  according  to  length  scale  associated  with  the 
discontinuity;  i.e.,  L.  The  three  parameters  have  the  follow¬ 
ing  form: 

1)  Nu  (Gr/4) =  -  H ' (0) (x  +  I)-*  (47) 

0 

2) '  NuJGr/4)'5*  =  -  H'  (0)  (x  +  1)  -  1 . 40e  *  (48) 

3)  Nu(Gr/4)  =  -  [H*  (0)  +  (Z  *  (0)  +  e*Z'  (0)5 

10  11 

+  fZ'  (0)  +  £*Z'  (0)  +  £*2Z’  (0) )  52  ]  (49) 

20  21  22 

The  heat  transfer  parameters  were  plotted  as  functions  of  x, 
i.e.,  for  suction  and  blowing  magnitudes,  £*,  of  .05  and 
.25  in  Figures  (26) -(29).  As  expected  in  the  region  from  the 
leading  edge,  x  =  -  L,  to  the  discontinuity  at  x  =  0  the  heat 
transfer  rates  for  cases  1)  and  3)  are  identical  and  in  case 
2)  the  rate  is  higher  by  a  fixed  amount  over  the  unperturbed 
free  convection  case.  At  the  discontinuity  the  curve  for  the 
discontinuous  case  experiences  an  abrupt  change  in  slope  and 
intersects  the  continuous  case  curve  at  approximately  ^  = 

0.22  for  £  *  =  0.05  and  at  ^  =  0.52  for  z*  =  0.25.  The  magni¬ 
tude  of  applied  suction  apparently  has  a  significant  effect 
on  the  location  of  the  crossover  point  where  the  discontinuous 
heat  transfer  rate  exceeds  that  for  the  continuous  case.  The 
behavior  of  the  discontinuous  heat  transfer  rate  as  £  increases 
does  not  appear  to  approach  an  asymptotic  suction  limit  charac¬ 
teristic  of  a  constant  boundary-layer  thickness.  However, 
since  the  inner  temperature  series  expansion  constructed  to 
approximate  behavior  near  the  surface  is  only  valid  in  the 
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Figure  26:  Comparison  of  Heat  Transfer  at  the  Surface  of  the 
Plate  with  Suction,  e*  *  0.05. 
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Figure  2  7:  Comparison  of  Heat  Transfer  at  the  Surface  of  the 
Plate  with  Suction,  e*  =  0.25. 
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Figure  28 


Comparison  of  Heat  Transfer  at  the  Surface  of  the 
Plate  with  Blowing,  e*  =  -0.05, 
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neighborhood  of  the  discontinuity,  i.e,,  small  values  of  jj, 
it  cannot  be  expected  to  accurately  predict  results  in  asymp- 
totic  limiting  case  as  -j-  becomes  large.  The  addition  of  higher 
terms  to  the  inner  expansion  would  serve  to  extend  the  range 
of  validity,  however,  since  the  Sparrow  and  Cess  analysis  of 
the  continuous  case  employed  only  first  order  terms,  it  is 
reasonable  to  expect  the  results  of  the  present  analysis  to 
provide  a  closer  approximation  to  the  physical  situation. 

This  essential  difference  between  the  two  studies  should  be 
kept  in  mind  when  comparing  results.  The  situation  with 
blowing  applied  as  depicted  in  Figures  (28)  and  (29)  shows  the 
expected  reduction  in  heat  transfer  rates  when  compared  to 
the  unperturbed  free  convection  in  both  the  continuous  and 
discontinuous  cases.  At  the  discontinuity,  however,  the  curve 
representing  discontinuous  blowing  exhibits  a  downward  shift 
in  slope  and  then  follows  a  path  which  is  basically  parallel 
to  the  free  convection  and  continuous  blowing  curves.  The 
crossover  that  occurs  in  the  case  of  a  blowing  magnitude  of 
e*  =  -0.05,  is  probably  an  anomaly  due  to  an  insufficient 
number  of  terms  in  both  the  continuous  and  discontinuous  series 
expansions.  For  a  blowing  parameter  of  £*  =  -0.25  all  three 
curves  appear  to  asymptotically  approach  the  same  slope.  A 
final  judgement  as  to  the  degree  of  accuracy  with  which  the 
present  analysis  models  the  heat  transfer  behavior  with  discon¬ 
tinuous  blowing  or  suction  applied  must  await  comparison  with 
experimental  data. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


An  overall  review  of  the  results  obtained  from  this  analysis 
leads  to  a  number  of  general  conclusions.  An  apparent  range  of 
validity  for  the  expansions  can  be  estimated  from  a  comparison 
of  the  inner  and  outer  profiles  presented  in  Figures  2-25  based 
on  the  degree  of  correlation  between  them.  In  terms  of  the 
streamwise  distance  from  the  discontinuity  a  qualitative  upper 
bound  of  j-  =  1.0  is  suggested.  In  terms  of  the  blowing  and 
suction  parameter  the  relationship  given  on  equation  (44)  sug¬ 
gests  an  upper  limit  of  |e*|=  1.0  which  equates  to  a  magnitude 
of  approximately  .05  umax  for  v  for  which  the  original  boundary 
layer  equations  should  remain  valid.  In  the  case  of  suction 
the  inner  profiles  appeared  to  model  the  expected  behavior  of 
the  location  and  magnitude  of  the  maximum  velocity  with  respect 
to  the  unperturbed  free  convection  profiles  better  than  did 
the  outer  profiles.  When  blowing  was  applied,  the  outer  pro¬ 
files  correlated  with  anticipated  behavior  in  this  respcet 
better  than  the  inner  profiles.  In  general  the  matching  of 
the  inner  and  outer  velocity  profiles  improves  as  the  stream- 
wise  distance  from  the  discontinuity  decreases  and  is  also 
qualitatively  b^-^er  in  the  blowing  cases.  The  analysis  also 
indicated  that  while  the  inner  series  provides  reasonably  good 
results  at  very  small  distances  from  the  plate,  in  an  overall 
integrated  sense  the  outer  profiles  exhibit  a  closer  correlation 


with  the  velocity  behavior  associated  with  free  convection. 

The  effects  of  suction  and  blowing  on  the  temperature 
distribution  both  appear  to  be  well  modeled  by  the  inner 
series  expansions.  Again  the  degree  of  matching  between  the 
inner  and  outer  profiles  significantly  improves  as  the  stream- 
wise  distance  from  the  discontinuity  decreases.  Also  the 
effects  of  the  discontinuity  on  temperature  behavior,  relative 
to  the  free  convection  case,  are  less  direct  than  on  velocity 
behavior.  This  difference  is  evidenced  in  the  manner  of  coup¬ 
ling  exhibited  by  the  ordinary  differential  equations  which 
define  the  inner  expansion  coefficient  functions  (see  equations 
(19)  and  (20)).  This  fact  in  addition  to  the  apparent  over¬ 
estimate  of  velocity  magnitudes  produced  by  the  outer  series; 
the  underestimate  by  the  inner  series;  and  the  heat  transfer 
behavior  predicted  for  small  suction  magnitudes  suggest  that 
additional  higher  order  terms  be  added  to  both  inner  and  outer 
expansions  to  achieve  a  better  approximation  of  the  physical 
situation.  In  general  the  results  obtained  for  the  effects 
of  discontinuous  suction  or  blowing  on  heat  transfer  at  the 
surface  of  the  plate  seemed  consistent  with  the  results  for 
unperturbed  free  convection  and  the  case  of  continuous  suction 
or  blowing  applied  over  the  entire  plate  length.  Comparison 
with  experimental  data  when  it  becomes  available  will  enable 
a  more  definitive  conclusion  to  be  made  regarding  the  validity 
of  this  solution  with  respect  to  heat  transfer  performance. 
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Also,  with  the  addition  of  higher  order  terms  to  the 
inner  series  it  would  be  possible  to  determine  if  the  heat 
transfer  results  predicted  by  the  present  analysis  tend  to 
approach  the  asymptotic  limits  for  large  values  of  £  es¬ 
tablished  by  Merkin  [9].  In  that  study  the  following 
asymptotic  limits  for  continuous  suction  and  blowing  at 
large  distance  from  the  leading  edge  were  reported,  given 
below  in  terms  of  the  parameters  used  in  this  analysis: 

1)  For  suction  Nu/fp-)^  -*■  4e*Pr  as  £  ■+•  °° 

2)  For  blowing  Nu/  (p)  4  -*•  0  a s  jj  -*■  00 

The  asymptotic  limit  in  the  case  of  suction  is  represented 
by  the  dashed  line  in  Figures  (26)  and  (27) . 

Although  an  exact  determination  of  the  inner-outer  series 
match  point  is  not  essential  when  the  primary  interest  of 
the  investigation  is  to  predict  behavior  at  the  surface  of 
the  plate,  a  more  rigorous  determination  of  that  point  would 
serve  to  enhance  the  overall  creditability  of  the  present 
method  of  solution.  Therefore,  another  recommendation  to 
improve  the  method  point  is  to  use  the  energy  balance 
integral  given  in  (42)  to  form  an  energy  difference,  i.e., 
from  a  difference  integral  by  substracting  the  right  side 
from  the  left  side.  A  numerical  minimization  procedure 
could  then  be  applied  to  this  difference  to  determine  the 
match  point  which  most  closely  satisfies  the  energy  balance. 
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APPENDIX  A 


CONSTRUCTION  OF  THE  FREE  CONVECTION  MATCHING  SERIES 
The  equations  which  define  the  behavior  of  the  free  con¬ 
vection  velocity  and  temperature  profiles  were  presented  in 
equations (9a)  and  (9b),  with  appropriate  boundary  conditions, 
in  terms  of  similarity  functions  F  and  H  which  characterize 
the  stream  function  and  temperature  distribution  respectively. 
In  the  region  near  the  surface  of  the  plate,  y*  =  0 ,  and  at 
the  location  of  the  discontinuity,  x  =  0,  these  functions  can 
be  approximated  by  suitable  series  expansions  in  powers  of 
y*  as  follows: 

F(y*)  =  Z  a.(y*)j  (A-l) 

j=o  J 

H(y*)  =  1  +  1  b.(y*)j  (A-2) 

j-1  3 

where  the  variable  y*  can  be  used  since  it  is  identically 
equal  to  the  free  convection  similarity  variable  at  x  =  0 . 

The  constant  coefficients  in  (A-l)  and  (A-2)  can  be 
determined  from  applying  the  following  previously  enumerated 
boundary  conditions: 

At  y*  =  0:  F ( 0 )  =  0 

F'  (0)  =  0 
H(0)  =  1 

(A-  3) 

As  y*  -*■  00 :  F '  (»)  =  0 

H(»)  =0 
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These  conditions  immediately  yield  the  result  that  a  =  a  =  0. 

0  1 

By  evaluating  the  series  representations  for  F"  and  H'  the 
constants  a  and  b  are  specified  as: 

2  i 

a  =  ?sF"(0)  and  b  =  H’(0) 

2  1 

where  both  of  these  initial  values  are  well  known  and  depend 
only  the  Prandtl  number. 

To  evaluate  the  higher  order  coefficients  the  form  for 
F ' ' '  must  be  obtained  from  equation  (9a): 

F' ' '  =  -  3FF"  +  2F ' 2  -  H 


which  when  evaluated  at  y*  =  0  yields  F '  '  ' (0)  =  -  1  and 

equating  this  value  to  the  series  form  of  F ’ ’ *  yields: 

1 

a3  =  '  S’ 

A  similar  technique  is  used  in  determining  b  from  equation 

2 

(9b)  and  the  following  value  is  obtained: 


b  =  0 
2 

The  remainder  of  the  coefficients  are  found  by  successively 
differentiating  equations  (9a)  and  (9b)  and  equating  them  to 
the  appropriate  series  representation  evaluated  at  y*  =  0. 
The  results  of  this  process  are  summarized  below: 


a  =0 
o 

a  =  0 

1 

a  =  *sF"(0) 

2 


a  = 

3 


1/6 


b  =  H’ (0) 
2 

b  =  0 
2 

b  =  0 

3 


(A-  4) 


b  =  - PrF" (0)  H  * (0) 

4  8 
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-H' CO) 
a  =  — 

4  24 

a  =  f F" (0)  ]  2 

5  120 

These  series  forms  for  F  and  H  will  be  used  to  satisfy  the 
matching  conditions  for  both  the  inner  and  outer  expansions 
later  in  the  course  of  the  analysis. 
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APPENDIX  B 


DETERMINATION  OF  THE  STREAM  FUNCTION  AND  TEMPERATURE  DISTRI¬ 
BUTION  SERIES  EXPANSION  PARAMETERS 

Recalling  that  the  definitions  of  the  stream  function  and 
temperature  profile  series  expansions  involved  certain  arbi¬ 
trary  constants,  the  methods  and  procedures  used  in  evaluating 
them  are  outlined  in  this  appendix. 


i 


*  =  e*(Z/<Pi)a  *  3 Cc [ Ckf £ Cn) ] 

(B-l) 

e  =  l  +  4>uCd[Ckgk(n)l 

(B-2) 

where 

a  K 

C  =  <j)  x  and  n  =  <P  yk* 

1  2 

(B-3) 

The  object  of  constructing  such  expansions  for  the  stream 
function  and  temperature  distribution  is  to  obtain  series  re¬ 
presentations  which  satisfy  the  momentum  and  energy  equations 
which  mathematically  model  the  physical  problem.  The  indivi¬ 
dual  components  of  those  governing  equations  (6)  and  (7),  i.e., 
the  partial  derivatives,  must  be  evaluated  in  terms  of 
expressions  (B-l)  through  (B-3).  Substituting  (B-l)  and  (B-2) 
into  the  momentum  and  energy  equations  yields: 

Momentum  Equation 

l-2b  2b-l 

4  4>  a  S2c  +  a  ( [a(c  +  k)  +  bKkf']Kkf’] 

12  3  K  k 
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,2b  c  +  2b 

[a(c  +  k)^kfk][Ck£^]}  -  a  **4.^  a  [£kf"] 


-  <t>  a  O3*  5  a  Kkfv'"3  +  1  +  4>  [Ckgv] 

12  3  K-  4  K. 


(B-4) 


Energy  Equation 

1-b  _ . j, b-1 

-  c+d+ -  ^  , 

4<fri  a  <t>J3<t>J  a  ([a(d  ♦  k)5  gk][CK£^] 

_b  b  +  £ 

-  [a(c  +  k)Ckfk3  [Ckg^]>  -  4e**i  [?kg£] 


_2b  2b_  +  d 


(B-S) 


Certain  conditions  which  must  be  satisfied  will  impose 
restrictions  on  the  allowable  values  for  the  <j>’s  and  constants 
a,  b,  c  and  d.  The  one  condition  which  must  be  satisfied  by 
the  equations  is  the  initial  velocity  and  temperature  matching 
at  a  distance  L  from  the  leading  edge  of  the  plate  where  the 
discontinuity  in  normal  velocity  begins.  At  this  location  the 
perturbed  velocity  and  temperature  profiles  must  match  the 
free  convection  profiles.  In  terms  of  the  stream  function  the 
velocity  is  given  by: 


*  34>  _  a  c,  "b  A  ry*^ 

u*  =  -r~*  *  <t>  <P  <P  C r— ) 


ak.  in.  a  _  xv  a. 

b  a  (?p)  b  fk]  (B'6) 


The  free  convection  velocity  profile  written  in  terms  of  u* 
can  be  expressed  as  follows: 


u*  =  F '  Ly  * ]  *  £  ja.  (y*)' 
j  =  2  J 


(B-  7) 
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Equating  the  leading  term  in  each  series  (B-6)  and  (B-7)  i.e.. 


when  j=2  and  k=0  we  obtain  the  following  relationship 

ca  ac  , 

c  -r —  -r  ' i 

2a  y*  =  4>  <j> 

2  1  2 


y3  n  J 


(B-8) 


For  proper  matching,  the  exponents  of  y*  must  be  equal;  therefore, 
ac 


’5~ 


-  1 


(B-9) 


Also,  since  it  is  desirable  to  have  expressions  which  are  in¬ 
dependent  of  e*  for  the  boundary  conditions  on  the  as  Ty-<», 

the  following  condition  will  be  imposed; 

_ca 

b  <b  =1  (B-10) 

1  2  3 

The  requirements  for  matching  the  temperature  profiles  yields: 
ad  _ka  _ka 

4>c4>  ^  (^)[4>k<}>  *”(£)  b  gk]  *  .^b-y^  (B- 11) 

i2  4  n  12  n  x  j-J-3 

Applying  similar  arguments  to  those  used  for  the  condition  of 
velocity  matching,  the  following  relationships  are  obtained 
from  (B-ll)  : 

-Jp  =  1  (B-12) 

and 

<$>  =1  ( B - 1 3 ) 

1  2  4 

For  additional  restrictions  on  the  exponent  constants  consider 
equation  (B-4)  in  which  the  highest  order  derivative  term  must 
appear  in  the  Oth  order  expansion  in  order  to  obtaine  meaningful 
results.  This  requirement  imposes  the  following  condition  on 
the  exponents  of  K- 


c  + 


3b 

a 


<  2c  + 


2b -1 
a 
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or  equivalently, 
ac  >  b  +  1 


(B-14) 


The  same  restriction  also  applies  to  the  energy  expansion, 
equation  (B-5)  and  the  same  relationship  is  obtained  (B-14). 
Since  fractional  exponents  on  £  are  not  permissable,  due  to 
the  inability  to  match  coefficients  of  like  of  powers  of  £  in 
both  expansions,  from  equation  ( B - 1 2 )  the  following  further 
information  is  obtained: 

d  =  1  (B- 15) 

and 

-  £  =  1  (B- 16) 

Mow  combining  (B-9)  and  (B-16)  ,  the  constant  c  is  specified 
as : 

c  =  2  (B-17) 

Using  this  information  and  the  relationship  given  in  (B-14) 
imposes  the  following  restrictions  on  the  possible  values  of 
a  and  b : 

a  >_  1/5  and  b  <_  -  1/3  ( B - 1 8 ) 

The  relationships  previously  obtained  for  the  <j>'s  can  be  sim¬ 
plified  by  substituting  the  known  values  for  the  exponents  in 
equations  (B-10)  and  (B-13)  from  which  are  obtained  the 
following : 


$2<|>2<j>  =  1 

(B-19) 

1  2  3 

d>  4*  =  1 

(B-20) 

1  2  M- 
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Since  equation  (B-5)  is  homogeneous  in  <j>  and  4>  appears  only 

4  4 

singly  in  equation  (B-4),  without  loss  of  generality  it  can 

be  set  equal  to  unity.  Setting  $  =  1  in  equations  (B-19)  and 

(b-20)  yields  the  following  further  relationships  for  the  <J>'s: 

<J>  <p  =1  (B-  21) 

1  2 

and 

4>  =  1  (B-22) 

The  foregoing  analysis  has  permitted  the  reformulation  of 
equations  (B-4)  and  (B-5)  in  terms  of  only  two  unspecified 
parameters  a  and  $  ,  and  those  equations  reduce  to  the  follow¬ 
ing  forms:  ,  , 

Momentum:  a{[k  +  1)  ?kf£ ]  Ukf£]  -  [  (k  +  2Kkfk)  • 

[Ckfy]  }  -  4e*CUkf£]  =  Ukf£’]  +  5  +52Ukgk]  (B  -  2  3) 

I  3.1 

Energy:  4^^  a{[(k  +  1)  Ckgkl  Ukf£ ]  *  [  (k  +  ZKkfk]  * 

-  4ennkg{]  =  p^rCCkg^3  ( B  -  24) 

Equation  (B-18)  allows  for  only  two  reasonable  choices  for 
the  value  of  a;  namely,  a  =  1/3  or  a  =  1/2.  Furthermore,  for 
each  of  these  two  values  <)>  can  be  defined  to  either  include 
e*  or  not;  thereby  posing  four  possible  cases  for  consideration. 
The  -cases  which  specify  a  =  1/2  do  not  satisfy  the  boundary 
condition  for  f  as  and  therefore,  are  unacceptable  for 

i 

the  present  problem.  The  case  for  which  e *  is  included  in  $ 
and  a  =  1/3  results  in  a  set  of  ordinary  differential  equations 


92 


which  include  e*  in  the  homogeneous  part  and  which  cannot  be 
transformed  into  general  equations  independent  of  e*  through 
the  use  of  universal  functions.  Therefore,  the  case  with 
a  =  1/3  and  z*  not  included  in  the  definition  of  <j>  was 

i 

selected  to  complete  the  solution  of  the  problem.  A  summary 
of  the  results  is  presented  below: 


1/3,  b  =  -  1/3,  c=  2 ,  d  =  1 

(B- 25) 

,1/3  ,  -  1/3 

(j)  ,  0  -  (j)  ,  4  =  <P  =  1 

4  2  4  3  4 

(B- 26) 

3  1/3  3  -  1/3 

(fx)  ,  n  -  (£x)  y*  -y/5 

CB- 27) 

j  e*C3  4  C 2  [ Ckfk ] 

(B- 28) 

1  +  C[Ckgk] 

( B -  29) 

The  momentum  and  energy  expansions  which  result  from  these 
definitions  were  presented  in  Chapter  II,  equations  (17) 
and  (18).  All  that  remains  is  to  specify  the  appropriate 
boundary  conditions  for  the  f^  and  gk . 

At  n  -  0:  u*  =  =  5Ukf£]  -  0,  (B-  30) 

V*  ='H  =  ~[e*  +  J  £*l[(k+2Kkfk]}  -  ±e\  (B- 31) 

and 

9  -  1  +  C [Ckgkl  =  1  (B- 32) 

From  these  original  boundary  conditions,  the  local  boundary 
conditions  for  fk  and  f£  are  determined  to  be: 

fk(°)  =  f'(0)  =  gk(0)  =  0,  for  all  k  (B-33) 

The  boundary  conditions  at  rr*-00  are  obtained  through  matching 
the  velocity  and  temperature  profiles  to  the  free  convection 
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profiles.  The  velocity  match  is  performed  by  equating  equations 
(B-6)  and  (B-7)  with  the  appropriate  values  substituted  for  the 
unspecified  parameters.  The  result  of  this  process  is  the 
following  relationship: 


j~2 


j  (y  *) 


j'x  -  n^)k+1f£] 


(B -  34) 


By  equating  like  powers  of  v*  and  substituting  the  values  for 
the  a^  obtained  in  Appendix  A  the  following  local  boundary  con¬ 
ditions  for  the  f£  as  rv*°°  are  determined: 
f ' 

— —  =  2a  =  F" (0)  as 

r,l  2 


f  ’ 

— i-  =  3a 

Ty  2  3 


1 

2 


as 


(B- 35) 


-i-  =  4a  -  -  as  n- 

_  3  4  O 


[F"(0)1 


J.  =  Sa  =  _ 

s  24 


as 


The  temperature  profile  matching  yields  the 


boundary  conditions  for  the  g^: 


following  local 


H'  (0)  as 


as  n-*00 


as  tv*-00 


PrF" (O)H' (0) 
- p-— v- 


as  n-*00 


(B- 36 ) 
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APPENDIX  C 


THIRD  ORDER  UNIVERSAL  FUNCTION  EQUATIONS 
The  third  order  ordinary  differential  equations  which 
result  from  the  expansions  given  in  (17)  and  (18)  are  pre¬ 
sented  below,  although  only  functions  of  second  order 
and  lower  were  employed  in  the  actual  analysis: 

f '  "  +  2f  f”  -  5 f ’  f'  +  5 f" f  =  Sf'f'  -  3f’f"  -  4f"f 

3  03  03  03  '12  12  12 


-  g"  +  2f  g’  -  4 f ' g  +  5g ’ f 

Pr  3  0*3  0S3  0  3 


g  f'  -  2g  f' 

0  3  13 


3g  f' 
2  1 


-  3f  ^  g ’ -  4f2gj  -  4e  *g '  (C-2) 

The  corresponding  forms  for  f  and  g^  in  terms  of  the  univer¬ 
sal  functions  and  the  resulting  differential  equations  are  as 
follows : 

f  =  Y  +  e  *Y  +  £  * 2  Y  s*3Y  (C-  5a) 

3  3  0  3  1  3  2  3  3 

and 

g=Z  +e*Z  +  e*2Z  +  e  * 3  Z  (C-3b) 

330  31  32  33 


The  third  order  momentum  equation  becomes: 

M  [Y  ]  =  5Y '  Y’  -  3Y  Y"  -  4Y"  Y  -  Z  ;  (C-4a) 

3  30  1020  1020  1020  10 

M  [Y  ]  =  5Yr  Y'  +  5Y’  Y’  -  3Y  Y”  -  3Y  Y" 

3  31  1120  1021  1021  1120 
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NL 


*7 


(C- 4b) 


-  4Y"  Y  -  4Y"  Y"  -  -  , 

1  3  2  1  1  1  2  0  1  1 


M  [Y  ] =  5Y’  Y '  +  5Y'  Y'  -  3Y  Y"  -  3Y  Y" 

3  32  1022  1121  1022  1121 


4Y"  Y  -  4Y"  Y  -  4Y"  ; 
1022  1121  21 


M  [Y  ]  =  5Y ’  Y ’  -  3Y  Y”  -  4Y"  Y  -  4Y" 

3  33  1122  1122  1122  22 


(C- 4c) 
(C-4d) 


where  the  linear  differential  operator  M  is  defined  by: 


M  «  Y,M  +  2f  Y"  -  5f 'y '  +  Sf"y 

3  3  0  3  0  3  0  3 


(C-4e) 


The  third  order  energy  equation  becomes: 


E 

[Z 

]  =  11 

Y* 

+ 

3Z 

Y  *  - 

3Y 

Z'  - 

4Y 

Z  ' 

y 

3 

3  0 

1 

0  2  0 

2 

0  1  0 

1  0 

2  0 

2 

0  1 

0 

E 

[Z 

]  =  11 

Y" 

+ 

11 

Y'  + 

3Z 

Y"  + 

3Z 

Y 

t 

3 

3  1 

1 

0  2  1 

1 

1  2  0 

2 

0  1  1 

2  1 

1  0 

3Y  Z* 

-  3Y 

Z' 

-  4Y 

Z' 

-  4Y 

Z' 

. 

4Z 

10  2  1 

1 1 

2  0 

2 

0  1  1 

2 

i  i 

0 

E 

[Z 

]  =  11 

Y’ 

+ 

11 

Y'  + 

3Z 

Y'  + 

3Z 

Y 

t 

3 

3  2 

1 

0  2  2 

1 

1  2  1 

2 

i  l  i 

2  2 

1  0 

_ 

3Y  Z' 

-  3Y 

Z' 

-  4Y 

Z' 

-  4Y 

Z' 

. 

4Z 

10  2  2 

1  1 

2  1 

2 

1  1  1 

2  2  1 

0 

E 

[Z 

II 

is) 

tNJ 

Y f 

+ 

3Z 

Y’  - 

3Y 

Z'  - 

4Y 

Z 

» 

3 

3  3 

l 

1  2  2 

2  2  11 

i 

1  2  2 

2  2 

i  i 

(C- 5a) 


2  o 


2  1 


(C-5d) 


where  the  linear  differential  operator  E  is  defined  by: 

E  =  i-Z"  +  2f  Z'  -  4f ' Z  +  5Z '  Y  -  Z  Y'  .  (C-5e) 

3  rr  3  03  03  003  003 

The  boundary  conditions  for  these  equations  are: 


at  n  =  0: 

Y 

3  0 

=  Y 

3  1 

-  Y 

3  2 

=  Y 

3  3 

=  0  ; 

Y* 

3  0 

=  Y' 

3  1 

=  Y' 

3  2 

=  Y' 

3  3 

-  o  ; 

(C-6) 

and 

Z 

3  0 

=  Z 

3  1 

=  Z 

3  2 

z 

3  3 

0; 
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and  as  ^  3  0  [F"  (0)  ] 2  .  y> 

n4  24  31 

and  h±  ■»  -PTF'(0)H»_(0I.  z  ^0; 

n4  8  3  1 

which  follow  from  the  coefficient 

given  in  equations  (B-33),  (B-3S) 


-*■  0  ;  Y '  -*■  0  ;  Y '  -  0 

3  2  3  3 

(C-7) 

Z  -*•  0;  Z  -►  0 

3  2  3  3 

function  boundary  conditions 
and  (B- 36) . 
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APPENDIX  D 


DETERMINATION  OF  THE  OUTER  SERIES  COEFFICIENTS 
The  process  of  determining  the  coefficient  functions  of 
the  outer  series  expansions  defined  in  equations  (29)  and 
(30)  involves  the  solution  of  the  ordinary  linear  differential 
equations  (32);  the  evaluation  of  the  resulting  constants  of 
integration;  and  the  solution  of  the  algebraic  equations  (33). 
The  procedure  is  greatly  simplified  if  the  outer  expansions 
for  4/  and  0  are  first  transformed  into  series  with  coefficients 
which  are  powers  series  in  y* ,  the  coefficients  of  which  are 
derived  from  the  asymptotic  forms  of  the  f^  and  gk-  This  tech¬ 
nique  will  facilitate  the  evaluation  of  the  constants  of  in¬ 
tegration  significantly  reduce  the  complexity  of  equations 
(39). 

The  asymptotic  forms  of  the  inner  series  coefficients  can 
be  written  as  polynomial  in  n  with  as  yet  undetermined 
coefficients  as  follows: 


f 

2 

3 

a  n2 

2 

(D-la) 

f 

3 

An3  +Bn2+Cn+D 

(D- lb) 

l 

i  iii 

f 

3 

An1*  +  B  n3  +  c  n2+  D  n  +  E 

(D-lc) 

2 

2  2  2  2 

2 

g 

=r 

b  n 

(D- Id) 

0 

i 

g 

= 

B*n  +  c* 

(D-le) 

l 

i  i 

g2 

* 

B*n  +  C*n  +  D 

2  2  2 

(D-lf) 
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which  must  satisfy  the  asymptotic  boundary  conditions  specified 
in  ( B - 35 )  and  (B-36).  The  constants  are  evaluated  by  substi¬ 
tuting  the  equations  (D-l)  into  the  original  differential 
equations  (19)  and  (20)  and  then  solving  the  algebraic  equa¬ 
tions  which  appear  as  coefficients  of  like  power  of  n.  The 
results  of  this  process  are  summarized  below: 


A  = 

l 

-1 

6 

A 

2 

_  _b  i 

24 

B*  = 

i 

0 

B  = 
i 

0 

B 

l 

=  0 

C*  = 
i 

biCi 

2a 

2 

(D-2) 

C  = 
i 

TBD 

C 

2 

.-Ci 

"  Ta~ 
2 

B*  = 
2 

0 

D  = 

l 

'I  c* 

3  e 

D 

2 

=  TBD 

c*  = 
2 

■b  iC ! 

8a2 

2 

E 

2 

C? 

4a 

D*  = 
2 

-b ,  D, 

6a2 

where  a 

2 

*  \  F"(0)  ,  b 

=  H» 

l 

(0) 

and  C 

and  D 

1  2 

to  be 

determined 

from  application 

of  the  matching  condition. 

This 

process 

will  be 

outlined 

at  the  end 

of 

this 

appendix. 

Now 

the  equations 

(D-l) 

with  the 

coefficients 

as  defined 

in  (D-2)  can  be  substituted  into  stream  function  (15a)  and 
temperature  distribution  (15b) .  Then  if  the  definition  of 
n  =  yVS  is  invoked  and  all  coefficients  of  like  powers  of  i 
grouped  together,  the  following  expansions  are  obtained  which 
are  valid  as  y*  becomes  large: 

*  *  [a2y*2-  i  y*3  -  ^  y*4]  +  [0]C 

+  [Cy*  -  *£2y*2K2  +  [D2y *] c3  (D-3) 
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and  9  =  [1  +  b  y*]  +  [0]5  + 


b  D  1 

.  i  2  r  3 


b  C  b  C 
*  x  -  *  1  v*  iz 

2*7  8a! 

2  •* 


CD-4) 


With  these  equations  in  hand  the  final  determination  of  the 

outer  series  coefficients  can  now  proceed. 

The  first  of  the  differential  equations  defined  in  (32a) 

is  logarithmic  and  its  solution  has  the  form: 

ip  =  K  ip  '  (D-5) 

1  1  o 

Since  ip  is  the  coefficient  of  the  linear  term  in  the  outer 

i 

expansion  and  that  coefficient  in  (D-3)  is  zero,  must  also 

equal  zero,  and  ^has  the  final  form: 

ip  =0  (D-6) 

i 

Using  this  equation  in  (D-2b)  reduces  it  to  logarithmic  form 


with  solution  for  ip  as  before: 

2 

tp  =  K  \p  ' 

2  2  0 


(D-7) 


where  K  is  evaluated  by  equating  (D-7)  with  the  appropriate 
2 

form  of  the  £2  coefficient  in  (D-3).  The  final  result  for 


ip  then  becomes: 
2  C 

■p  =  -d-  p  • 


(D-8) 


The  last  of  the  differential  equations  (32c)  can  be  expressed 
in  a  simpler  form  using  (D-6);  however,  its  solution  requires 
a  different  technique.  The  reduced  form  for  (32c)  is: 


i  \p  ’\p  ’  -  i  tp  >p  "  =  0  +  ip  ' 

£  0  3  ^  3  0  0 


CD-9) 
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Using  the  free  convection  stream  function  equation  (9a)  p " ' 

0 

can  be  written  as: 


r 


ip "  +  Zip  '  2  —  9 

0  9  0  0 


(D- 10) 


and  by  suitable  algebraic  manipulations  in  combining  (D-9) 
and  (D-10)  the  following  differential  form  is  obtained: 


2  (D- 11) 

o  /  \r  o  ' 

Straight  forward  integration  yields  the  following  general 


solution  for  \p 


ip  =  6\p 
3  0 


2y*p  '  +  K  ip  1 


(D-12) 


where  is  obtained  as  before  with  the  aid  of  equation  (D-3). 

Thus  the  particular  solution  has  the  following  form: 

D  3 

P  =  6\p  -  2y*p  ’  +  ^  ip  ’  (D-13) 

30  o  ^a  o 

2 

This  completes  the  determination  of  the  outer  expansions 

for  ip  and  9  with  the  exception  of  the  constants  C  and  D  . 

1  2 

Since  these  constants  appeared  originally  in  the  asmyptotic 

forms  of  f  and  f  and  are  determined  from  the  matching  of  the 
1  2 

inner  and  outer  series,  the  results  of  equating  the  universal 
function  forms  to  the  asmyptotic  forms  will  yield  expressions 
for  C  and  D  .  The  following  equations  result  from  this 


process : 
For  f 

j 

and  for  f 


-in3  *  C  n 
O  0  10 


t  £  *  =  Y  +  e  *Y 
0  10  11 


(D- 14) 


■ki.  n* 
24  o 


C_l 

4a 


nz  +  D  n  + 

0  2  0  4a 
2  2 


Y  +  e  *Y  +  e*2Y 
20  21  22 


(D-1S) 
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(D - 16) 


! 


Solving  for  C  from  (D-14)  yields  the 

C  -  C  +  e*C 
,  10  11 


following : 


where 

C 

i  o 


and 


=  ( — LL  +  *) 

li  n  3nJ 


A  similar  procedure  yields  the  following  form  for  D  obtained 

2 

from  (D-15)  : 


D  =  D  +e*D  +  e  *2 
2  20  21  22 


D 


CD-17) 


where 


D 

2  0 


b 

i 

24 


n 3+ 


C 

dr 

2 


c2 

n  '  TaV  ’ 

2 


D 

2  1 


and 


Y 


C  C 
io  m  . 

2a  n  }’ 
2 


D  =  (_il  -  _LL  ) 

22  n  4a  n 

2 

where  the  value  n  was  taken  as  nen(^  the  chosen  end  point  of 

numerical  integration.  This  completes  the  derivation  of  the 

outer  series  expansions,  except  to  note  that  the  constants 

C  ,  C  ,  D  ,  D  ,  and  D  are  evluated  numerically  from  the 
10  11  20  21  22 

tabulated  universal  functions. 
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APPENDIX  E 


INTEGRAL  ENERGY  BALANCE  ANALYSIS 
In  only  heat  transfer  problem  an  overall  thermodynamic 
energy  balance  must  be  satisfied.  For  the  case  of  the  present 
problem  the  integral  energy  balance  takes  the  following  form: 

/°°c  P  (T  -  T  )  udy  =  /X+Lg"  (x)dx  (E-l) 

op  0  0 

The  left  hand  side  reduces  to  the  following  form  with  non- 

dimensional  integrand: 

„  ^  Gr  1/4  - 

4CpPv(rw  -  TJ(JI)  /  u*6dy*  (E-2) 

Since  there  are  two  distinct  regions  over  which  this  integral 

is  defined,  it  must  be  split  into  two  integrals  separated  by 

the  point,  y^gg,  where  the  inner  and  outer  expansions  should 

match.  Substituting  the  appropriate  forms  for  u*  and  9  as 

defined  by  the  inner  and  outer  series  the  integral  (E-2)  can 

be  written  as  the  sum  of  the  following  two  integrals: 

Inner  region  y*  =  0  to  y*  =  y, 

4  cDyp  (T  -T  )  ( )1/4  /Y{f*(nK  +  [f'(n)g  (n)  +  f'(n)H2 

♦  [f'(n)g  (n)  +  g  (n)f’(n)  +  f’(n)]C3}  dy*  (E-3a) 

0  1  0  1  2 

where  Y  is  the  matching  point, 

V* 

n  =  ,  and  dy*  =  £dn  which  transforms  (D-3a)  into  the 

following : 

<f;«2  *  tf,,8,  *  f;i«!  * 
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(E- 3b) 


[f'g  +  g  £'  +  f']C3}dn 
0  1  0  1  2 


where  y  depends  on  x. 

For  the  outer  region,  y*  =  y  to  y*  *  the  following  integral 
applies : 

i  c 

4cppv(Tw-T)(^)%  /°°*F’H  +  jL  (F’H'  ♦  F'HK2 

2 


T  D  1  1  D 

+  Wlk  -  ^  +  J  F’H  +  6 


-  2y*) F"H] 5 3 )dy*  (E-4) 


The  right  hand  integral  in  (E-l)  also  must  be  divided  into 
two  distinct  regions,  the  first  extending  from  the  leading 
edge  to  the  discontinuity  and  the  second  from  the  discontinuity 
to  any  arbitrary  location  above  it.  The  following  two  inte¬ 
grals  are  obained: 

(E-S) 


/Lq"dx  +  q"  dx 
on  o  nP 


where  qp  is  the  heat  flux  in  the  undisturbed  free  convection 
region  of  the  plate  and  qp  is  the  heat  flux  above  the  discon¬ 
tinuity,  both  evaluated  at  the  surface  of  the  plate.  Upon 
making  suitable  substitutions  to  non-dimensionalize  the  in¬ 
tegrands,  the  following  result  is  obtained  from  (D-4) : 

-  4k(Tw-Tj(^-)!5  [i  H '  (0)  +  |  g'C3  +  \  g’  C4 


+  r  g'?5l  ( E - 6 ) 

D  2 

where  Gr  is  evaluated  at  L  and  the  g£  are  evaluated  at  the 
surface  of  the  plate,  n  =  0. 

After  completing  the  indicated  integrations  in  (E-3a)  and 
(E-3b)  and  equating  coefficients  of  like  powers  of  C  the 
following  relationships  are  obtained: 
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H'(0)  -  H'(Y)  =  3PrF(Y)H(y)  (E-7) 

C^F'  (y)H(y)  -  [FM (0 )  ]  2Y  2  CH-8) 

i-F"  (O)H'  (0)  (^— -)  +  f  (£)  -  F(y)H(y)  -  H’(y) 

j  ^21?,  3Pr 

D  1  1 
-  jjj  F’(y)H(y)  +  jYF'(y)H(y)  =  -  H’(0)  (E-9) 

2 

The  relationships  can  be  identically  satisfied  only  at  the 
surface  of  the  plate,  n  =  0 .  From  this  fact  it  can  be  concluded 
that  the  outer  profiles  are  the  better  approximation  of  global 
behavior  and  will  satisfy  the  energy  balance  in  an  integrated 
sense,  pointwise  local  behavior  is  required  near  the  surface 
of  the  plate  the  inner  expansions  should  be  used,  although 
they  cannot  satisfy  the  integral  energy  balance. 
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